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Abstract 



We consider a massless scalar field obeying Dirichlet boundary conditions 
on the walls of a two-dimensional L x b rectangular box, divided by a movable 
partition (piston) into two compartments of dimensions a x b and (L — a) x b. 
We compute the Casimir force on the piston in the limit L — > oo. Regardless 
of the value of a/b, the piston is attracted to the nearest end of the box. 
Asymptotic expressions for the Casimir force on the piston are derived for 
a <^.b and a>i». 
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I. INTRODUCTION 



In 1948 Casimir predicted a remarkable macroscopic quantum effect: two conducting and 
neutral parallel plates should attract each other due to the disturbance of the vacuum of the 
electromagnetic field caused by their presence [1] (for a general review on the Casimir effect, 
see Ref. [2]). Inspired by that result, a few years later Casimir suggested that the zero-point 
pressure of the electromagnetic field might yield the stresses postulated by Poincare in order 
to explain the stability of the electron [3]. Boyer, however, showed that the Casimir force 
for a conducting spherical shell is repulsive [4], thus invalidating Casimir's model for the 
electron. 

Boyer's result brought attention to the fact that the attractive or repulsive character of 
the Casimir force depends on the geometry of the configuration. This has been investigated 
in detail for fields (scalar or electromagnetic) confined in a rf-dimensional rectangular box 
[5-11]. Let us consider, for instance, a massless scalar field subject to Dirichlet boundary 
conditions at the walls of the two-dimensional box < x < a, < y < b. The vacuum 
energy is formally given by (h — c — 1) 



One can perform the summation using analytic regularization (AR); the result is (see Ap- 
pendix A) 



where Z 2 is an Epstein zeta function [12]. An analysis of (2) shows that the sign of the 
Casimir tension T = —dE 0jAR /dA (where A = ab is the area of the box) depends on the 
ratio b/a: it is positive if 1 < b/a < 2.74 and negative if b/a > 2.74 [2]. 

There are, however, at least two reasons for which one should be suspicious of the use 
of Eq. (2) as the basis for such an analysis. First, it does not take into account the con- 
tribution to the vacuum energy from the region outside the box, which, in principle, also 
depends on its dimensions. (This problem was discussed recently in [13], but the solution 
proposed there is incomplete.) Second, its finiteness is an artifact of the AR scheme: more 
than just regularizing integrals or sums, it also does a certain amount of renormalization 
by automatically subtracting power-law divergences (in this respect, AR is similar to di- 
mensional regularization [14]). This is precisely what happens here. If one regularizes the 
sum over modes in Eq. (1) with a smooth cutoff function and performs the sum using the 
Abel- Plana formula, one obtains (see Appendix B) 



with Ci(A) ~ A 3 and C 2 (A) ~ A 2 as A — > oo. (We have discarded terms that vanish in that 
limit.) 

The difference between -Eo, cutoff and -Eo,ar would be harmless if the first two terms on 
the r.h.s. of Eq. (3) could be absorbed into counterterms. Let us forget for a moment the 



E (a,b) = u jk , 
j,k=i 




(1) 





E , cutoff (a, b) = Ci(A) ab + C 2 (A) (a + b) + E 0tAR (a, b), 



(3) 
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problem of neglecting the exterior modes, and examine this question. The first term has the 
form eQab, where eo is the energy density of the vacuum in the absence of the box. It can be 
cancelled by a "cosmological constant" counterterm, a constant added to the Hamiltonian 
density in order to make the vacuum energy in free space equal to zero. The problem lies in 
the second term: being proportional to the perimeter of the box, it may be interpreted as 
(part of) the self-energy of its walls. Such a term cannot be eliminated by a renormalization 
of the parameters of the theory [15,16]. (This problem also occurs in the parallel plates 
configuration. In that case, however, it can be ignored if one is interested only in the force 
between the plates, for their self-energies do not depend on the distance between them. 
In the present case, the dismissal of the self-energy of the box walls could be justified if 
perimeter-preserving deformations are the only ones allowed.) 

In this work we shall examine a slightly different system in which both problems can be 
ignored. Instead of the box discussed above, we shall consider a box of dimensions L x b 
divided by a movable partition, or piston, into two compartments, A and B, of dimensions 
a x b and (L — a) x b, respectively (see Fig. 1). If one is interested — as we are — in 
computing the Casimir force on the piston, then the contribution to the vacuum energy 
from the region outside the box can be ignored, as it is not affected by the position of the 
piston. In addition, as will be shown below, the divergent terms in the Casimir energy are 
naturally eliminated when one computes the force on the piston. We shall compute this 
force in the limit L — > oo and show that it pulls the piston to the nearest end of the box 
regardless of the value of the ratio a/b. We shall also derive asymptotic expressions for the 
force for a <C b and a^> b. 



II. CASIMIR FORCE 

The total energy of the vacuum for the system described in the previous paragraph (and 
depicted in Fig. 1) can be written as the sum of three terms: 

E = E£ + E* + E™\ (4) 

Using the cutoff regularization discussed in Appendix B, the first two terms are given by 
E$ = E 0: cutoff (a, b) and E$ = E 0:CUtoS (L - a,b) [see Eq. (3)], so Eq. (4) becomes 

E = E 0t AR (a, b) + E , AR (L -a,b) + d(A) Lb + C 2 (A) (L + 2b) + £ out . (5) 

The Casimir force on the piston is given by —dE /da. Since the last three terms on the 
r.h.s. of Eq. (5) do not depend on the position of the piston, we obtain the following result 
for the Casimir force on it: 

F = ~ [E , AR (a, b) + E 0yAR (L - a, b)} . (6) 

As anticipated, although the total vacuum energy contains divergent terms and a term (E™ 1 ) 
that one does not know how to compute, the Casimir force on the piston is finite and can 
be computed exactly. 
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The result one obtains for the force inserting (2) into Eq. (6) is not very illuminating, 
so, before we actually compute F, let us derive an alternative expression for Z? ,ar(«) b). In 
order to do that, it is convenient to define the auxiliary function 



S(m,a;s):=n-^T(-) £ 



m\ ■ 



7T / 



+ 



n 



-s/2 



(Re(s) > 1). 



(7) 



Its analytic continuation to the complex s-plane (with simple poles at s = 1, —1, —3, . . .) is 
given by [7] 



S(m, a; s) = 



am 



n (l-s)/2 



S - 1 



+ 4E 



if(i_ s )/ 2 (2nma) 
(nmfl)( 1_s )/ 2 



(8) 



where K v (z) is the modified Bessel function. Eqs. (7) and (8) allows us to reexpress the 
Epstein zeta function that appears in Eq. (2) as [17] 



oo 

Z 2 (a,6;3)= £ ' (j V + k 2 b 2 ) 



-3/2 



j,k=— oo 
oo oo 



= E' E (iV + ^ 2 )- 3/2 + E'(^ 2 ) 

j=—oo k=— oo 



-3/2 



fe=— oo 



27r3/2 V* <tf ■ 1 /I. * a. 2C ^ 

_g5(7rja,l/6;3) + -^- 

2tt 2 16tt ~ k T ^ ( , a\ 2C(3) 
+ TI2- E 7^ + 



3a 2 6 a& 2 , ? 

Inserting this result into Eq. (2) yields 

vr C(3)a 1 



6 3 



(9) 



00 k 

3,k=l •> 



(10) 



Inserting Eq. (10) and the corresponding expression for _E ,ar(£ — a, b) into Eq. (6) and 
taking the limit L — > oo we obtain the following expression for the Casimir force on the 
piston: 



67' 



(11) 



where K[(x) = dK\[x)jdx. Since i^i(x) is a monotonic decreasing function of x, it follows 
from Eq. (11) that F < for all (positive) values of a/b; in other words, the piston is 
attracted to the nearest end of the cavity. 

It is easy to obtain an asymptotic expression for F valid for a ^> b: since K\(x) ~ 



7r/2x exp(— x) for large x, one may retain only the term with j — k — 1 in Eq. (11), thus 
obtaining 



F „ JL ( a6 3)-V2 eX p 



2iia 



(a > 6). 



(12) 
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This result has the same form as the asymptotic expression of the Casimir force between two 
plates in one-dimension in the case of a scalar field with mass m = n/b [7]. This fact has a 
simple physical interpretation: when o>5 the system becomes quasi-one-dimensional, with 
the field acquiring an effective mass equal to the energy gap A = n/b due to the confinement 
in the transverse direction. 

In order to obtain an approximation to F valid for a < 6, we note that ^o, ar(o, b) = 
Eo,AR,(b,a), so that we can replace Eq. (10) by 

E ^-Wa-lB-ipJ K ^)- <13) 

If, on the other hand, we still express E ^(L — a,b) in Eq. (6) according to Eq. (10), we 
obtain an alternative expression for the force on the piston (in the limit L — > oo): 

The last term in Eq. (14) is exponentially supressed when a <C b, so in this case we have 

„ C(3)6 , vr C(3) 



8yra 3 48a 2 16tt6 2 



(a < 6). (15) 



If one divides both sides of Eq. (15) by b, the first term on its r.h.s. correctly reproduces 
the Casimir tension between two infinite parallel lines a distance a apart [7]. The other two 
terms are subdominant for a <^b, and yield finite size corrections to that result. 



III. CONCLUSION 

We argued in this work that the knowledge of the vacuum energy inside a rectangular 
cavity is not enough for one to calculate the Casimir force on its faces. Two ingredients are 
missing in such a calculation: the knowledge of the contribution to the vacuum energy from 
the region outside the cavity, and the proper handling of divergent terms in the regularized 
expression of the vacuum energy. We then considered a slightly different type of cavity, 
namely, a rectangular box divided by a piston into two rectangular compartments. In this 
case, if one is interested only in the Casimir force on the piston, those ingredients can be 
neglected. In addition, the force-on-the-piston problem has two attractive features: (i) it is a 
simple generalization of the single-cavity problem, for which results are already available in 
the literature [5-11], and (ii) from the experimental point of view, it is simpler to construct a 
cavity with a piston than a variable-size rectangular cavity. Results for the electromagnetic 
field in a three-dimensional rectangular cavity with a piston will be presented elsewhere. 
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APPENDIX A 



Let us evaluate the divergent sum over modes in Eq. (1) using analytic regularization. 
We start with the function 



S(a,b;s) := ^ ]T 
z j,k=l 



-s/2 



(Al) 



which is defined for Re(s) > 2. As we shall see, its analytic continuation to the complex 
s-plane is well-defined at s = —1, so we can define the analytically regularized Casimir 
energy as E 0:AR (a,b) = £(a,b; -1). 

In order to obtain the analytic continuation of S(a, b; s) it is convenient to rewrite Eq. 
(Al) as [17] 



£(a,b;s) = | ]T' 

j=— oo h = — oo 



oo oo 



a) + ( b , 



-s/2 




-s/2 



^ [a) k tAK 



(A2) 



where Z p (a\, . . . ,a p ;s) and ((s) denote the Epstein and Riemann zeta functions, respectively. 
Applying the reflection formulae [7] 



r(0vr-n( S ) = r(^)vr^)/ 2 c(i- s ), 



(A3) 



Oi • • • a p r ( — ) 7T 



. . . , a p - s) = T 7r(-P)/ 2 Z p (l/ ai , . . . , 1/opjp 



(A4) 



to Eq. (A2) and taking s = — 1 we obtain Eq. (2). 



APPENDIX B 



In this Appendix we derive Eq. (3) using the Abel-Plana summation formula [18], 
1 r°° r°° F(e + it)-F(s-it) 



X; ^(n) = - F(0) + / Fit) dt + i lim / 
£r Q 2 Jo e^o+ Jo 



e 2nt _ 1 



dt. 



(Bl) 



F(z) is an analytic function in the right half-plane, going to zero sufficiently fast as \z\ — > oo, 
|axg(*)| < tt/2. 

In order to apply the Abel-Plana formula to the series (1), we have to introduce a smooth 
cutoff function D^(z,w): 
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n oo / -2 pX 1 / 2 



a' 6 . 



(B2) 



The function D A (z, w) must satisfy the following conditions in the region Re(z), Re(w) > 0: 
(i) it is analytic in both variables; (ii) it is real for z and w real; (iii) it vanishes sufficiently 
fast for \z\, \w\ — > oo (so that the regularized series is absolutely convergent); (iv) it is 
symmetric, i.e., D A (z,w) = D A (w,z), and (v) liniA^oo D A (z, w) = 1. An example of such a 
function is given by D A (z,w) = d A (z) d A (w), with d A (z) — [1 + [z + 1) 2 /A 2 ]~ 2 . 

Applying formula (Bl) to the series Sj in (B2), we can rewrite each of them as a sum of 
three terms, namely, 



(B3) 



S 



(2) 



(B4) 



: (3) 



i lim 

e^0+ Jo 



du 



2nu 



' j 2 u 2 .Y /2 n (j e + iu\ 



(B5) 



Applying (Bl) to £• yields 



i=i 7o 



V D ( V 

2a A \a' 



1 x 
+ - 



v dv 



a Jo e 



2wv 



IV 



D\ I , ) + c.c. 



(B6) 



Changing the variable of integration in the first integral to t = v/a and taking the limit 
A — > oo in the second one we obtain 



oo /*oo 



24a 



(B7) 



Similarly, application of (Bl) to X^j 5'| 2 ' ) yields 



VS 2) = -- dttD A (0,t) + ab du dv(u 2 + v 2 ) l/2 D A (u,v) + I A (a,b), (B8) 
~^ 2 jo Jo Jo 



j 

where 



/*oo /"C 

I A (a,b) = iab lim / dw / 

e^o+ Jo JO 



dv 



o2nav 



[u 2 — v 2 + ie) 1/2 D A {e + iv,u) - c.c. . (B9) 



Taking the limit A — > oo in the integral above we obtain 
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fOO pOO H^j 

lim Ua, b) = -lab / du / — — (v 2 - u 2 ) 1 ' 

a^oo AV ' Jo Ju e 2nav - 1 v ; 

= _2a6 [°° [ V du (v 2 - u 2 ) 1 ' 2 

Jo e 27Tav -1 Jo v ' 

nab r°° v 2 dv 
Jo 



2 Jo e 2nav - 1 
C(3)6 



8n 2 a 2 



(BIO) 



Finally, let us take the limit A — > oo in Eq. (B5). Changing the variable of integration -u to 
£ = au/jb, we obtain 



a^oo J a 2 ii e 2 ^ bt / a - 1 

00 /'OO 

7 1 ^ 1 



2l oo 



(7 2 

U fc=l 



= "— E T f 2^' -V (Bll) 

Collecting all pieces together we obtain Eq. (3), with E AR (a,b) in the form given by Eq. 
(13) and 

C 1 (A) = - du dv(u 2 + v 2 ) 1/2 D A (u,v), (B12) 
2 Jo jo 

roc 

C 2 (A) = -n dttD A (t,0). (B13) 
Jo 
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FIGURES 



b 


A 


B 




a 


L — a 



FIG. 1. Two-dimensional L xb rectangular box. A movable partition (a piston) divides it into 
two compartments, A and B, of dimensions a x b and (L — a) x b, respectively. We shall assume 
that a, b <C L — > oo. 
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